Non-Markovian Stochastic Liouville equation and its Markovian representation. 
Extensions of the continuous time random walk approach. 
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Some specific features and extensions of the continuous time random walk (CTRW) approach 
are analyzed in detail within the Markovian representation (MR) and CTRW-based non-Markovian 
stochastic Liouville equation (SLE). In the MR CTRW processes are represented by multidimen- 
sional Markovian ones. In this representation the probability distribution function (PDF) W(t) 
of fluctuation renewals is associated with that of reoccurrences in a certain jump state of some 
Markovian controlling process. Within the MR the non-Markovian SLE, which describes the effect 
of CTRW-like noise on relaxation of dynamic and stochastic systems, is generalized to take into 
account the influence of relaxing systems on statistical properties of noise. The generalized non- 
Markovian SLE is applied to study two modifications of the CTRW approach. One of them considers 
the cascaded CTRWs in which the controlling process is actually CTRW-like one controlled by an- 
other CTRW process, controlled in turn by the third one, etc. Within the MR simple expression for 
the PDF W(t) of total controlling process is obtained in terms of Markovian variants of controlling 
PDFs in the cascade. The expression is shown to be especially simple and instructive in the case 
of anomalous processes determined by long time tailed W(t). The cascaded CTRWs can model the 
effect of complexity of a system on relaxation kinetics (in glasses, fractals, branching media, ultra- 
metric structures, etc.). Another CTRW-modification describes the kinetics of processes governed 
by fluctuating W(t). Within the MR the problem is analyzed in a general form without restrictive 
assumptions on correlations of PDFs of consecutive renewals. The analysis shows that W(t) can 
strongly affect the kinetics of the process. Possible manifestations of this effect are discussed. 

PACS numbers: PACS numbers: 05.40.Fb, 02.50.-r, 76.20.+q 



I. INTRODUCTION 

Relaxation and transport phenomena in condensed me- 
dia are of great importance for understanding of various 
processes in physics, chemistry, biology etc. [U, 3- 
For description of these phenomena a large number of 
approaches have been proposed in literature 0, H, 0, @ . 
The most popular are Markovian approaches based on 
Langevin and Fokker-Plank 0, 0] equations assuming 
negligibly short memory in processes under study. Great 
advantage of these approaches consists in possibility of 
relatively simple treatment of relaxation in dynamical 
and (Markovian) stochastic systems induced by interac- 
tions whose fluctuations are Markovian processes (MPs). 
This type of relaxation, which in what follows for brevity 
is called MP affected MPs, is known to be described by 
the stochastic Liouville equation (SLE) [5|. 

Recent investigations show, however, that Markovian 
theories do not properly describe the kinetics of relax- 
ations and fluctuations of various observables in many 
non-equilibrium processes [f| 0] • The non-Markovian re- 
laxation kinetics in these processes results form strong 
memory effects. Moreover, in a number of systems the 
kinetics is anomalously long time tailed which implies 
long time memory in them. The anomaly manifests it- 
self in some interes ting phenomena such as non-er godicity 
0,113, aging slow relaxation 0, M H Ef, 

etc. 

Anomalous specific features of relaxation in highly 
non-equilibrium systems attract close attentions of in- 



creasing number of scientists 0, 0, 13- Analysis of 
these features can be made within different models and 
approximations. One of the most popular and, probably, 
adequate is continuous time random walk (CTRW) ap- 
proach in which the process is represented by jump-like 
fluctuations (or relaxation) treated as a sequence of re- 
newals [3, 0, HH characterized by the probability density 
function (PDF) W(t) of waiting times of renewals [this 
PDF is often denoted as ip(t) 

Last years different variants have been analyzed and 
applied to treating a large number of processes 0, EH • 
In some recent studies special attention have been paid 
to above-mentioned anomalously long time tailed relax- 
ation 0, E3 : within the CTRW approach such behavior 
is assumed to result from that of the PDF W(t) ~ l/t 1+a 
with a < 1 0, Qj|. The applications of the CTRW ap- 
proach allowed one to interpret a lot of experimental re- 
sults 0,E3]- It is shown, in particular, that in many cases 
the specific properties of fluctuations and relaxation can 
quite properly be described within the anomalous CTRW 
approach (with anomalous behavior of the PDF W(t)). 

The CTRW approach is fairly popular last years. It is 
successfully applied both to general investigations and to 
interpretation of experimental results. It is worth point- 
ing out, however, that in majority of most interesting 
applications the validity of this approach is not quite ev- 
ident. Typically, it is justified with the use of intuitive ar- 
guments concerning statistical properties of systems un- 
der study. In addition, some recent investigations show 
that, despite great flexibility, the conventional CTRW 
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approach proves to be not quite accurate in description 
of the kinetics of a number of non-Markovian processes 
[l2l [l9| . In such a case any further extensions of this 
approach are of course desirable. 

Certain modifications of the conventional CTRW ap- 
proach have already been discussed in literature. In some 
CTRW- variants the non-homogeneity of the process, con- 
sisting in the dependence of the jump-like fluctuations 
on the fluctuation number, have been taken into ac- 
count [2(| [2l[ ■ In other variants the modification of sim- 
ple CTRW time sequences of renewals is proposed [HJ . 
These modifications are very interesting and essentially 
clarify specific features of CTRW-like processes. 

In this work we will propose and analyze some ex- 
tensions of the CTRW app roach within the Markovian 
representation 0, [23|, [24]]. This method is based on 
idea that under fairly general assumptions a wide va- 
riety of non-Markovian CTRW processes can be treated 
as MPs with fluctuating parameters, whose fluctuations 
are described by other MPs called hereafter controlling 
processes. In so doing these CTRW processes are un- 
ambiguously represented by some multidimensional MPs 
[20I ] and the statistics of renewals is associated with that 
of reoccurrences in certain transition state (or the state 
of onset of fluctuation jumps) during the controlling MPs 
in some auxiliary spaces. In particular, the PDF W(t) is 
expressed in terms of the characteristic functions of the 
controlling processes. 

The Markovian representation is shown to be very use- 
ful for the analysis of non-Markovian CTRW-like mod- 
els describing different physical processes [2(J HE HH- 
In addition it permits rigorous derivation of the non- 
Markovian stochastic Liouville equation (SLE) treating 
relaxation in dynamical and Markovian stochastic sys- 
tems induced by CTRW- t ype fluctuating interactions 
with thermal bath [2(| l23l l24j. In what follows, for 
brevity, Markovian systems in which relaxation is in- 
duced by CTRW-fluctuating interactions will be called 
CTRW affected MPs. 

In this work the compact formulation of the Markovian 
representation is proposed which simplifies and general- 
izes derivation of the CTRW-based non-Markovian SLE 
thus allowing for extension of the SLE to take into ac- 
count back effect of fluctuating system on effective sta- 
tistical properties of fluctuations. This formulation espe- 
cially clearly demonstrates that CTRW processes can be 
considered as MP affected MPs (in above-proposed brief 
terminology) and basic equations of the CTRW approach 
are nothing else but the SLEs in a reduced form. 

The obtained generalized non-Markovian SLE is ap- 
plied to the analysis of validity of CTRW approaches. 
Two important extensions of the CTRW approach are 
proposed and discussed in detail: 

The first modification describes the effect of cascaded 
controlling processes in which the process controlling re- 
newals is assumed to be CTRW-like one controlled by the 
second CTRW process, which in turn is controlled by the 
third CTRW processes, etc. In this cascaded CTRW con- 



trol model the compact representation for the PDF W(t) 
(more correctly for the Laplace transform of this func- 
tion) is obtained in terms of PDFs of controlling pro- 
cesses Wj(t) at all cascade steps j, found assuming these 
processes to be Markovian (i.e. neglecting control) . This 
model is very suitable for the analysis of non-Markovian 
relaxation kinetics in structured and disordered systems 
especially in the case of anomalous long time tailed be- 
havior of PDFs Wj(t). 

The second modification treats CTRW-like processes 
governed by fluctuating PDFs W(t). This type of pro- 
cesses can be considered as an extension of conven- 
tional CTRW-approach in which fluctuations of W(t) 
are assumed to result from the additional effect of non- 
equilibrium medium and are modeled by the dependence 
of the system on the Markovian stochastic variable of 
special type. In this model and within the Markovian 
representation the description of these processes reduces 
to solving the non-Markovian SLE. The analysis of the 
model shows that fluctuations of the PDF Wit) can 
strongly change the kinetics of CTRW relaxation. Espe- 
cially significant effect is expected in the case of anoma- 
lous long time tailed processes. 

II. GENERAL FORMULATION 

We consider relaxation processes in a dynamical or 
stochastic Markovian system induced by fluctuating in- 
teraction with the classical thermal bath. The Markovian 
evolution of the system is assumed to be governed by the 
fluctuating Liouville operator L(t). Our general formula- 
tion is applicable both to quantum and classical systems. 
In particular, in the dynamic systems the only difference 
of these two cases is in the form of the operators L: for a 
quantum system L = i [H, . . . ] , (here [H , . . . ] is the com- 
mutator with the Hamiltonian H of the system), while 
for a classical system L = {H, . . . } (with {H, . . . } be- 
ing the classical Poisson brackets). In what follows, for 
definiteness and brevity, we will concentrate on general 
results as applied to classical systems. 

Fluctuations of L(t) are assumed to result from the 
dependence on the fluctuating bath coordinate x: L{t) = 
At(t)) whose changes are modeled by stochastic jumps 
between states \v) = |x„) with different L = L v in the 
multidimensional space {x}. 

Hereafter we will use " bra-ket" notation for the states 
in {x}-space suitable for treating relaxation phenomena 
determined by not self-adjoint evolution operators. 

In the model under study the system evolution is de- 
scribed by the (PDF) p(t) (or density matrix for quantum 
systems) which satisfies the linear equation 

p = —Lp with L = L x = y |x I/ )ii / (x !/ |. (2.1) 

For simplicity of presentation, L is considered to be diag- 
onal in {a;}-space though the formulas obtained are valid 
in the case non-diagonal L as well. 
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The operator solution of eq. (|2.1[) is represented as 



p(t) = U(t)p with U(t) = f 



(2.2) 



where T is the time ordering operator and po = p(t = 0) 
is the initial condition. 

Experimentally measured observables are usually de- 
scribed by the evolution operator U(t) averaged over 
stochastic fluctuations of x(i) which is expressed in terms 
of the conditional evolution operator G(x,x'\t) as 

(U(t)) x = V 6(x,xo|t)«7i(xo) = (x e |G| Xl ), (2.3) 

' — 'x.xo 

where <Ji(x) = |xi) is the initial PDF of the system in 
{x}-space normalized by the condition ^ x er,;(x) = 1. In 
eq. (|2.3p we have also introduced the (adjoined) equilib- 
rium state vector which in bra-ket notation is represented 
as [see below eq. (|3 . 8[) ] : (x e | = £)„{x„|. 

In accordance with eq. (|2.3|) the problem reduces to 
evaluating the operator G(x, Xj|t). Unfortunately this 
can be done relatively easily only for very few models of 
stochastic x(£)-fluctuations. In this work we will analyze 
some of those models based on the CTRW-approach. 



fluctuations as a sequence of sudden, jump-like changes 
of L. In the simplest variant of the CTRW (more com- 
plicated variants are discussed below) the onset of any 
particular change of number is described by the prob- 
ability Pp-i(t) (in {xj-space) not to have any change 
during time t and its derivative W@-i(t) = — P^_i(t), 
i.e. the PDF for times of waiting for the change. These 
functions are independent of 8 for 8 > 1 and for 8 — 1 
depend on the problem considered [H, 0] : 



W p > 1 (t) = W n (t), W {t) = Wi(t), 



(3.4) 



and Pj(t) — L dr W 3 ;(r), (j — n,i),. In what follows we 
will mainly discuss the non-stationary CTRW variant in 
which W > o {t) = W n {t). 

The Laplace transforms of Wj (t) andPj(i), (j — n,i)), 

are related by simple equation Pj(e) = [1 — Wj(e)]/e with 



W 3 {e) = [l + %{e)]-\ (j = n,i). 



(3.5) 



In eq. (|3.5p $>j(e) is the important auxiliary function [2C 



Markovian representation 



III. MODELS OF FLUCTUATIONS 
A. Markovian models 

The Markovian approach is based on the assumption 
that L(£)-fluctuations are described by the stochastic 
MPs in {x}-space and L(i)-evolution is described by the 
PDF er( X ,xo|i) satisfying equation 



cr = -Co- with er(x, x |0) = <5 x .x , 



(3.1) 



where C = C x is some linear operator in {a;}-space. The 
principal simplification of the problem results from the 
fact that in the Markovian approach (|3.1[) . i.e. in the 
case of MP affected MPs, G(x,x |t) obeys the SLE: 

G=-(L + C)G, so that G(t) = e - {L+t)t . (3.2) 

which yields for (U) = f™dt (U(t)) cxp(-ei) : 

0) = <G) X = (xeKe + L + r)- 1 ^}. (3.3) 

In the form (|3.2[) the SLE is valid for any dependence 
L x on the coordinate x though, in general, it is still very 
complicated for analysis either numerical or analytical. 

Significant simplification can be gained within some 
special models, for example, in the CTRW approach. 



B. CTRW-based models 

Non-Markovian L(£)-fluctuations can conveniently be 
described by the CTRW approach 0, 0, It treats 



The results of our earlier studies [20|, |23j show that im- 
portant non-Markovian generalizations of the SLE (|3.2j) 
can be obtained by assuming the operator C(t) to be a 
stochastic function of time. 

In this work we will analyze quite natural generaliza- 
tion of the Markovian SLE (|3.2[) based on the assumption 
that fluctuations of C{t) are Markovian. In other words 
C(t) = >C z ( t ) is assumed to be a function of the Markovian 
stochastic (controlling) variable z(t) which, in general, is 
a vector, so that generalizing the solution (|3.2p of the 
Markovian SLE one can write 



G(t) = f {exp[- Jodr(L x(T) + £,«)]}• 



(3.6) 



In the major part of the analysis we will assume that 
both {x}- and {zj-spaces are discrete and use bra/ket 
notation with Greek and Latin characters for states in 
these two spaces, for example, Ix^) and \zj), respectively 
(although, in case of need the continuous variants of {x}- 
and {z}-spaces will also be considered). 

We start our analysis with the simple Markovian 
stochastic migration process in {x}-space, which is de- 
scribed by the Kolmogov-Feller jump matrix 



C = k 



<t) -P X k% t) With P x = T,Uyt,lP^\ X »)( X P 



(3.7) 



in which = k d [z(t)] and = k n [z(t)) are the ma- 
trices (diagonal in {ir}-space) of z-dependent, i.e. fluctu- 
ating in time, jump rates and p vlJL are the probabilities of 
jumps {a;}-space normalized by the relation ^2 u Pvfi = 1- 
For simplicity, we assume that P x is independent of z, 
though most general results obtained in this section are 
valid in the case of z-dependent P x as well (see below) . 



4 



The matrix C(t) describes relaxation in {xj-space to 
the equilibrium state 

|Xe) = E^PSIXv), With (Xe| = EvW, ( 3 ' 8 ) 

for which (1 - P x )fc|x e ) = and (x e |(l - P x )k = 0. This 
state is assumed to be independent of z. Noteworthy is 
that even in the absence of the equilibrium state |x e ), for 
example when C describes diffusive migration in infinite 
space, the adjoint vector (x e | defined in eq. (|3.8|) still 
exists and satisfies the relation (x e |(l — P x )k = which 
means nothing else but the conservation of population in 
the process of migration in {x}-space. 

Within the considered Markovian approximation for 
£ z (t) -fluctuations the evolution of the system in control- 
ling {z}-space is governed by the PDF ip(z,z \t) satisfy- 
ing equation 

ip = -kip with y>(z|0) = \<fi(z)) = | Zi ), (3.9) 

in which A = A z is some linear operator describing relax- 
ation in {z}-space and <fi(z) = \zi) is the initial state of 
the controlling process (in this {z}-space) in the par-ket 
representation. Similar to the operator C in {x}-space 
for our further applications it is convenient to introduce 
the equilibrium states of the operator A 

K) = Ej#l»i) and = Eifal- (3-io) 

As in the case of {ir}-space considered above, the vector 
(z e | is still meaningful even in the absence of the equilib- 
rium state of the operator A. 

The operator A can, in general, depend on x. More- 
over, A can be a non-diagonal matrix in {x}-space 
though z for simplicity of further discussion, we will as- 
sume A to be diagonal in this space. 

In our analysis, in case of need we will use the simple 
and flexible continuous Smoluchowski model for control- 
ling process in {z}-space in which 

A = -^^-"'e^^Vj/^'e'WVj, (3.11) 

with z = |z|, u(z) is the effective potential, and V z = 
d/dz, is the radial part of the operator describing diffu- 
sion in n z -dimensional {z}-space with the diffusion coeffi- 
cient D z . This model is quite sufficient for our qualitative 
and semiquantitative analysis. The correspondence be- 
tween continuous model implying smooth functions D z 
and u(z) and its discrete variant considered above is for- 
mulated as follows: z <-> b\zj), where b is the spacing in 
the discrete model. 

Similarly to the simplest model discussed in Sec. III. A 
one can easily see that in this approximation the evolu- 
tion operator G r (r, r,-|i), which determines {U(t)) obeys 
the Markovian SLE in the extended space {r} includ- 
ing the dynamical subspace {x} and additional stochastic 
subspace {z}: {r} = {x ® z}: 



so that we get for G and the Laplace transform U = (G) x - 

G = (z e |G r |zj) and (U) = ((& + £ + A) _1 ) r , (3.13) 

where fl = e + L. In the form (|3 . 1 3|) the SLE is valid 
for any dependence of coupling C z on the coordinate z 
though, in general, it is still very complicated for analysis 
either numerical or analytical. 

Significant simplification can be gained within some 
special models for jump-rate matrix k z [see eq. (|2.3[) ]. 
One of the most well known is the CTRW-model dis- 
cussed below, in Sec. IVB. 

2. Localized jump rates. N on- Markovian SLE 

Very important results can be obtained in the model of 
interaction C z highly localized in {z}-space, i.e. highly 
localized jump rates k z and k z . The most interesting 
for our further discussion is the simplest variant of this 
model, which in the case of discrete {z}-space is written 
as 

k d , = fc P z ", K - hP? with Pi 1 = \z 3 ){z t \ (3.14) 

and j = t, n. Here kg is the matrix diagonal in {x}-space, 
while |z„) and |z t ) are the states {z}-space. The model 
(|3.14[) approximates £(i)-fluctuations by the set of jumps 
(blips) with the relaxation effect of each jump described 
by /co-matrix. As to properties of jumps statistics, they 
are determined by the controlling operator A and the 
localized matrix P™ 1 which implies that jumps occur as 
long as the system appears at the state |z t ) and each 
jump results in the transition |z t ) — > |z„) in {z}-space. 

Evidently, for the case of continuum {z}-space in eq. 
(I3.14|) [with A, for example, of type of that given by eq. 
(13. lip ] the term P™' should be replaced by the corre- 
sponding (5-function type one. 

Within the model (|3 . 14[) general formula (|3.13[) can be 
represented in a more suitable CTRW-like form: 

G(e) = (z e |G(e)[l-P,fc G(e)]- 1 |z 4 ) 

= P 4 (e)+P„(e)[l-W„(e)]- 1 W l (e), (3.15) 

where 

G(e) = {e + L + ki + k)- 1 , (3.16) 

is the Green's function, that essentially determines two 

matrices: the effective PDF-matrix Wj(e) of sudden 
fluctuation jumps (analogs of renewals) and the matrix 

Fj (e) of probabilities not to observe any fluctuation until 
time t, playing important role in CTRW-based theories 
0) S E3] [ see ec l- (|3.4p ]. These two matrices are ex- 
pressed in terms of the Green's function of the controlling 
process in {z}-space: 



G r = -(L + C + A)G r with G r (r,ri|0) = £ rr4 , (3.12) 



W J -(e) = P c 5? i (e), (j = i,n), 



(3.17) 
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in which 

fi*(e) = kiztlGie)^) = [1 + ^(e)]" 1 , (3.18) 

where 

%(e) = [g tt (e)ff tJ - 1 (e) - 1] + [feoS^)]" 1 (3.19) 

with 

g t] = (z t |g( e )| Zj ) and g{e) = (e + L + Ay 1 , (3.20) 
and 

Fj(e) = (z e |G(e)| Zj ) = (e + Ly^l - ^(e)], (3.21) 

In eq. (|3 . 2 1 1) {z e \ is equilibrium ket-vector in {z}-space 
(with (z e |A = 0) dehned by formula similar to eq. (|3.8p 
but with x-vectors replaced by those in {z}-space. 

Expression (|3.15[) . which in what follows will be called 
the generalized non-Mar kovian SLE (in resolvent form), 
looks similar to the simpler one derived earlier within the 
Markovian representation [2(J]. Moreover, for L = it 
reduces to that obtained in refs. [2(| [H[ in which the 
PDF matrices 

W,(e) = W Qj (e) = P x W 3 (e), (3.22) 
Fj(e) = P 0j -(e) = e- 1 [l-^(e)], (j=i,n), (3.23) 
where 

Wj(e)= di(z e |GoWI^)e" e ' = (z e |Go(e)|z,) (3.24) 
Jo 

is the PDF matrix in which (%(t) = e~^ fc ^ +A ' ) * and 

Go(e) = (e + k% + A) -1 , are the Green's functions de- 
scribing the stochastic process in {z}-space. Equations 
(I3.22p and (]3 . 23[) can be considered as a generalized vari- 
ant of relation (|3.5[) . 

Despite the similarity of obtained formulas with those 
known in the CTRW theory 0, 0], there are, however, 
some important specific features of the general expression 
(13. 15)) as opposed to the conventional CTRW formulas 
(PP ]) - (|531) and rf3"m (|3"2~3l . 

1) The matrix 

P 3 .(e) = ^- 1 (e)[l-fc (z t |(r!(e) + fco + A)- 1 |z,)] (3.25) 

in which f2(e) = e + L, can be non-diagonal, unlike simi- 
lar matrix in the conventional CTRW approach (see be- 
low). Moreover, in general, elements TPj ik (t) of the ma- 
trix Pj(t) do not satisfy the evident relation of type of 
Vj ik (0) = 1 which would allow one to interpret these ele- 
ments as probabilities. In addition, Pj(t) can be complex 
values, for example, in quantum processed for which the 
elements of L are complex. 

2) The relation between W,-(e) and Pj(e) implied by 
eqs. (|3.17l) and (|3.2ip does not coincide with that be- 
tween the PDF-matrix of waiting times and the matrix of 



probabilities known in the CTRW theory. Strictly speak- 
ing, the conventional relation [see eqs. (|3.22j) and (|3.23[) ] 
is reproduced only in the evident case L = 0. 

3) Equation (|3 . 1 5|) is of the form of CTRW expression 

however with matrices Wj(e) and P/(e) strongly modi- 
fied by the dynamic operator L [in contrast to the non- 
Markovian SLE appealing to the conventional PDF ma- 
trices Woj(i) and Pq,-(i) |c| (see eq. QEZB) )]. 

The above-mentioned effects of dynamic evolution can 
result in significant change of the time dependencies 
Wq/(t) and Poj(t) thus leading to the strong change of 
the kinetics of relaxation processes under study. 



3. CTRW-based non-Markovian SLE. 

The simplified variant of formula (|3 . 1 5[) can be ob- 
tained in the special case of matrix Wj(e) diagonal in 
{a;}-space, which is realized when the matrices A, L, and 
ko commute with each other: [A, L] = [L, ko] = 0. In this 

case one can obtain the representation for Wj (e) in terms 
of the Laplace transform of the conventional PDF-matrix 

of fluctuation blips Wq,- [see eqs. fl3~2"2"]) ) and (|3~2"4"1) )]: 

W j (e)=W Qi (n(e)) = P x dtW 3 {t)e- h ^\ (3.26) 
Jo 

where f2(e) = e + L. 

The simplified non-Markovian SLE (|3. 26[> . is neverthe- 
less more general than the variant of this equation in 

which Woj(e) independent of ^-coordinate (i.e. Wqj(e) 
proportional to the unity matrix in {x}-space) prj| . 

Noteworthy is that, in general, if A, L, and k z do not 
commute with each other, the representation (|3.26p is not 
valid and one should use the original expression (|3.15[) . 



IV. EXTENSIONS OF CTRW MODELS 

The proposed Markovian representation, based on de- 
scription of CTRWs as MP affected MPs, enables one to 
significantly extend the CTRW approach and the non- 
Markovian SLE. It reduces the treatment of system evo- 
lution to averaging exponential functional over Marko- 
vian fluctuations and offers the expression of CTRW 
equations in terms of multidimensional Markovian SLE. 

There are some straightforward extensions which do 
not need detailed analysis. For example, the evident vari- 
ant of extension is the model of several highly localized 
jump states in {z}-space (z ti | and |z ni ) (1 < i < i z , i z > 
1) in which k z = YsjkfPnt with P$ = |z„)(z t |. In 
this model the general formula (|3 . 13f> is also simplified 
by reducing the problem to solving the linear equation 
for the matrix G with elements Gjj(e) = (z ti \G(e)\z nj ). 
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Unfortunately, in this model the expression for G is fairly 
cumbersome and not quite suitable for applications. 

In this section we will mainly discuss less evident and 
more general extensions which allow for describing ef- 
fects of Markovian and non-Markovian fluctuating inter- 
actions on kinetics of relaxation in CTRW-like systems, 
i.e. (in our brief terminology) MP and CTRW affected 
CTRWs. In the Markovian representation the problem is 
still reduced to the analysis of MP affected MPs though 
in multidimensional space. 

Because of large number of parameters in this type of 
multidimensional MPs it is practically senseless to dis- 
cuss the problem in general. For this reason, here we 
will restrict ourselves to consideration of two modifica- 
tions and extensions of the CTRW approach interesting 
for theoretical analysis and for applications. 



A. CTRW-control. Coupled Markovian processes. 

1. General results 

So far in our analysis we have discussed the Marko- 
vian model for C(t) fluctuations which suggests that the 
controlling stochastic process z(t), responsible for fluctu- 
ations of the jump rates [see eq. (|3.7|) ] , is Markovian. 

Here we will consider the extension of the Markovian 
model based on the assumption that the controlling pro- 
cess z[t) is represented as a sequence (cascade) of control- 
ling Markovian processes in the multidimensional space 
{Zi} = {zi, Z2, • ■ ■ , zn} with evolution in each space {z q } 
being controlled by the process in {z 9 +i}-space whose 
mechanism is similar to that responsible for CTRW-type 
motion in {xj-space as described above in Sec. III.B. 
This mechanism implies that stochastic jumps in {z q }- 
space happen as long as the system appears in the tran- 
sition state \z q+1 ) (in {z g+ i}-space). As compared to the 
mechanism presented in Sec III.B, however, here, for sim- 
plicity, we will assume that jumps are not accompanied 
by the change of state in {z g+ i}-space, i.e. the final state 
z™ +1 ) = |z* +1 ). For the same reason we will also assume 
that initial state \z q ) of evolution in {z g }-subspace coin- 
cides with the transition state: |z*) = |z*) (1 < q < N). 

The above-formulated model of cascaded controlling 
processes is described by the following jump operators in 
spaces {Z q } = {z q , z 2 , ■ ■ ■ , z N } (q > 1) 

EN „ „ , 

AiPtK, where A, = (l-V^ko, (4.1) 
i—q 

with 

Pf=|z*)(z?|. (4.2) 

Here V Zi is the matrix of distribution functions of jump 
lengths (in {zi}-space) and fco i is the matrix of jump rates 
diagonal in {zi}-space. 



Similar to the model discussed in Sec. III.B, the 
control of motion in {xj-space will be described by z- 
dependence of the operator C: 

L = (l- P x )*oPf . (4.3) 

Noteworthy is that in accordance with obtained results 
the model (14. 1 [1 - (14. 3[) can be considered as a generalized 
variant of CTRW-type models for the controlling process. 

Recall that the problem under study reduces to eval- 
uating the evolution operator <G(t) defined in eq. (|3.6[) . 
In the proposed model of cascaded controlling processes 
this operator can be found in analytical form. 

Formula for G(t) can be obtained with the use of the 
general expression (|3.15|) 

G(e) = e- 1 [l-w(e)][l-P,u;(e)]- 1 (4.4) 
= e-^ie^+C^ 1 }- 1 (4.5) 

in which C x = (1 — P x )ko, 

w(e) = [1 + ^(e)]- 1 with $(e) = k- 1 ^ 1 (e+t 2 ) (4.6) 

and 

Me) = (zH^eK)- 1 = (z*|( e + AO^Izl)- 1 . (4.7) 

The expressions (14.41) - (|4. T[) relates the evolution opera- 
tor G(e) of the total system in the combined space {x®z} 
to that <h(e) = (e + Ai) _1 in the subspace {Zi}. 

The procedure presented above allows one to express 
4>i (£+£2) in terms of the evolution operator 32(e) in the 
reduced subspace {Z 2 }. The expressions similar to eqs. 
(|4.4[) - (|4.7|) . as applied to the operator ^"(e), yield 

0i(e) = (z 2 |[02(e + 4)+A 2 ]- 1 |z 2 )- 1 (4.8) 

with 

Ue) =gZ{z)- 1 (4-9) 

and 

92(c) = (4M*)\*l) = (4l(e + A 2 )- 1 |z 2 ). (4.10) 

By continuing the proposed procedure one gets the ex- 
pression for the function $(e) in terms of <7"(e) (with 
q > 1), and therefore in terms of functions 

^(6) = (z;|( e + A 9 )- 1 | z ;)- 1 , (4.11) 

which describe memory effects in the system as a result 
of Markovian motion in {z g }-subspace, i.e. without con- 
trolling interaction with other subspaces {z^ with i > q: 

He) = kQ 1 MMM- ■■)))&■ (4-12) 
2. Examples 

To illustrate the obtained results we will consider two 
limiting examples of two-state exponential and anoma- 
lously slow inverse power type controlling processes. 
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a. Two-state controlling processes. The exponential 
two-model for controlling processes corresponds to the 
simple expression for the evolution operator 

Me) ~ (z^ + A,)- 1 ^)- 1 

= e q /e- Kq /(e + Z q ), (4.13) 

where e q , K q , and t; q are some constant parameters. Ob- 
viously, formula (|4. 12|) with 4> q (e) presented in eq. (|4.13[) 
gives the expression for $(e) of type of continued fraction 
which predicts multiexponential behavior of the PDF ma- 
trix w{t). 

b. Anomalously slow controlling processes. Another 
very important model of controlling processes describes 
anomalously slow inverse-power type behavior of the 
PDF wi(t). This model is realized by taking a weak 
fractional-power type dependence of 4> q (e) 

4> q {e) = w q {e/w q ) a ", a q < 1, (4.14) 

where w q = Cg^o, and ( q ~ 1. Substitution of eq. (|4.14p 
into formula (|4. 12|) yields 

njv 
a q < 1 (4.15) 
q—l 

and Rq = Y[ q =i( w ql w q-i) lq ^ where j q = ]\ q i=1 on. 

Note that in the simplest case of identical rates w q : 
w q = w, (q > 1), when Ro = 1 and wn ~ w, formula 
(I4.15P is represented in a simple form $(e) ks (e/w) ao . 

B. Fluctuating CTRW-jumps 

1 . General formulas 

Another interesting problem in the CTRW theory, 
which can thoroughly be analyzed with the use of the pro- 
posed Markovian representation, concerns the kinetics of 
CTRW processes governed by stochastically fluctuating 
PDF matrices W/s(t) of jumps (renewals) or correspond- 
ing matrices wp(t) (here /3 is the number of jump). 

The fact is that in conventional CTRW theories the 
PDF Wp(t) and the probability Pp(t) are assumed to 
be the same for all renewals j3 > 1 except the first one 
(}3 = 0) and fixed functions of time [eq. (|3.4[) ]. In this 
section we will discuss the extension of CTRW processes 
in which the fluctuation jump kinetics is controlled PDF 
and probability matrices whose functional form fluctu- 
ates leading, in particular, to the difference of Wp(t) for 
different (3, 

The main difficulty in modeling fluctuating PDF ma- 
trices consists in necessity to take into account the nor- 
malization condition f£°dtWp(t) — Wp(e = 0) = 1, 
which ensures the population conservation in the process. 

The Markovian representation, discussed in Sec. Ill, 
allows us to make the above problem tractable. Within 



this representation the fluctuations of PDF functions 
Wf)(t) are assumed to be determined by the Markovian 
controlling process which is governed by the fluctuating 
operator A. The representation offers quite natural and 
fairly simple way of description of fluctuating A(i) by sug- 
gesting this operator to depend on the parameter y(t) 
which undergoes stochastic Markovian fluctuations. In 
this model the validity of the normalization relation for 
wp{t) can be ensured by taking proper form of the op- 
erator A(t). It is clear from eq. (|3.9p that this relation 
will be fulfilled if during evolution in {z}-space, described 
by fluctuating A(i), the population is conserved, i.e. eq. 
(I3.9|) can be represented in the form <p = — (V z -J z (<p, t)), 
where J z (ip, t) is the fluctuating flux in {z}-space. For in- 
stance, the population is certainly conserved in processes 
governed the Smoluchowski operator p. lip . More clearly 
this fact will be illustrated below with some examples. 

In general in the Markovian model proposed above, the 
PDF cr(y\t) representing j/(i)-fluctuation process satisfies 
equation of type of (|3.9p : 

a = -A y a with er(y,0) = cr, ; (y), (4.16) 

where A y is the operator responsible for the evolution in 
{y}-space. 

Similar to the case of L(t) fluctuations caused by 
Markovian fluctuating z(£)-parameter, the consideration 
of the effect of y (^-fluctuations reduces to analyzing the 
SLE for the evolution operator G r (£) in the combined 
space {r} = {x ® z ® y}: 

G r = -(£ + A z +A y )G r with G r (r,r 4 |0) = 6 rrz . (4.17) 
so that 

G = (y e \(z e \(e + £ + A y + A;)- 1 ^)^), (4.18) 

where |yj) = <xs(y) is the initial y(t)-state in bra/ket 
notation and 

(U) = ((e + C + Ay+KX 1 ),, (4.19) 

In eq. (|4.17p both the jump operator £ and the control- 
ling operator A z can, in principle, depend on variable y 
as a parameter: £ = £ v z and A z = A|. 

Just these dependencies of £ and A z on y model the 
effect of y(t) evolution [i.e. fluctuations of PDF functions 
Wj(t)] on the kinetics of processes under study. In gen- 
eral, both dependencies can strongly contribute to the 
effect, however in this work we will concentrate on the 
analysis of the effect of £|-dependence only, assuming 
that A z is independent of y. 

In the Markovian model of ilij (i)-fluctuations the prob- 
lem reduces to the analysis of the SLE (|4.17p which is sim- 
ilar to eq. (|3.12p considered above in detail. In general, 
eq. (I4.17P can be solved only numerically. In this work we 
will restrict ourselves to discussing some most important 
specific features of manifestation of ^ (^-fluctuations in 
the kinetics of the process within the simple analytically 
solvable variant of the proposed model. 
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2. Localized rate of fluctuation jumps 

The variant of the Markovian model which allows for 
analytical analysis of the problem can be considered as an 
extension of the approximation of highly localized jump 
rates discussed above. Here we assume high localization 
of fluctuation jump rates in both spaces {z} and {y}, i.e. 
in the combined space {z®y}. In this limit the expression 
for /^-dependence is given by 

CI = (l-PjkoP* with P*l = |y t z t )(y t z t |, (4.20) 

where |ytz t ) = |y t )|z t ) is the jump state in {z(g>y}-space 
which, for simplicity, is suggested to remain unchanged 
after jumps. 

In the model (|4.20|) y(t)-fluctuations are represented 
as a set of jumps (blips). The statistics of jumps is char- 
acterized by the PDF w y (t) of waiting times of jumps 
whose Laplace transform can be written as: 

€> y (e) = [l + $ y (e)}- 1 , (4.21) 

where 

*«(e) = [*d(yt|(e + A,)- 1 |y t )]- 1 . (4.22) 

It is easily seen that from mathematical point of view 
the proposed model is close that formulated in Sec. III. 
The difference consists only in a larger dimensionality of 
the space of the controlling process ({z(g>y}-space instead 
of {z}-one in the model considered in Sec. Ill) and in the 
corresponding difference in the form of the operator A 
which describes the controlling process (A z + A y instead 
of A,). 

This similarity enables us to use the results obtained 
above. It follows from these results that the model (|4.20[) 
predicts the CTRW expression for the evolution operator 

G of type of eq. (|3 . 1 5[) : 

G(e) = Pi(e) + P t (e)[l - W t (e)]- 1 W,(e), (4.23) 
where 

W,(e) = P x ti fj (e) and P,(e) = e" 1 ^ - G fj (e)], (4.24) 
(j = i,t), and 

4(e) = [l + $ / ,(e)]- 1 ! (4.25) 

however, with matrices 

H ( £ ) = [*o(y t z t |(e + A y + k z y l \y t ^)]-\ (4.26) 

whose behavior essentially depends on specific features 
of stochastic evolution of y(i). 

This expression makes it possible to easily analyze 
some specific features of the effect of y(i)-fluctuations. 
First of all, it is of special interest to discuss two natural 
limits: slow and fast y(i)-fluctuations. 



a. The limit of slow y(t)- and wj(t) -fluctuation. In 
the limit of slow fluctuations y(t), when the characteristic 
time of evolution in {y}-space, r y , is larger than that of 
the PDFs Wj(t), t w : r y 3> t w , one can neglect the term 
A y in eqs. (|4. 17(1 - fl4. 19|) thus arriving at the expression for 

G and (U) of type of those obtained in CTRW approach 
(see Sec. Ill) but averaged over y-distribution |<7j(y)) = 
|yi). In the short range model considered in this limit 

one gets the expression (|4.23[) - (I4.26[) for G(e) with: 

*/,(<0 « Ikaiytly^iztKe + K)- 1 ^)]- 1 . (4.27) 

b. The limit of fast y(t)- and uij(t) -fluctuations. In 
the opposite limit of fast fluctuations y(t), correspond- 
ing to T y < the SLE flUHD still yields the CTRW-like 
formula for the evolution operator G. In this formula, 

however, the matrices W,-(e) and Pj(e) are expressed in 
terms of the PDF matrix Wj(t) obtained by means of 
relations (|3.17|) - (|3.21|) which are evaluated with transi- 
tion matrices fcf(y) ana fc™(y) [see eq. (|3. 141) ] averaged 
over the equilibrium y-distribution cr e (y) = |ye)- I n the 
fast fluctuation limit the short range model results in the 

expression (|433>(B~2"6"j) for G(e) with: 

« [ko(yt\ye)(z t \(e + A z )- l \^)]-\ (4.28) 

where |y e ) is the equilibrium state in {y}-space. 

It is seen from formulas (|4.27p and (|4.28p that in the 
considered simple model of localized jump rates the ki- 
netics of the process is similar in both limits. The differ- 
ence is only in the characteristic parameters. However, 
in general, in the case of delocalized jump rates, i.e. delo- 
calized dependence of C v z on y, the kinetics in these two 
limits is, of course, different though the asymptotic be- 
havior at long times is expected to be similar. Below we 
will discuss the effect of delocalization of this dependence 
(see next Section). 

In addition, strictly speaking, in the above analysis of 
two limits we have assumed that y(t)-fluctuations are 
stationary which implies the existence of the equilib- 
rium state y e ) of the operator A y . It is of special in- 
terest to study the manifestation of non-stationary z(t)- 
fluctuations which can be fairly strong especially in the 
case anomalously slow fluctuations [tJ • 

3. Anomalous y(t) -fluctuations 

Here we will analyze in detail the effect of anomalously 
slowly fluctuating coordinate y(t). In the Markovian 
representation this type of fluctuations can be described 
with the use of the free diffusion model for ?/(i)-evolution 
discussed in Sec. III.B1. In this model the operator A y 
is given by 

A, = -y^VyiDyy^Vy), (4.29) 



9 



where V y = d/dy and n y is the {y}-space dimensional- 
ity (it is assumed that n y < 2). 

The model (|4.29p predicts the anomalous long-time- 
tailed dependence of the characteristic PDF w y (t) [see 
formulas P~2"Tj) and !|I7Z2"j) ]: 

«;,(*) - l/t 1+n y/ 2 . (4.30) 

These anomalous fluctuations of PDFs wp(t) of CTRW- 
jumps very strongly affect the kinetics of CTRW process. 
The effect can clearly be revealed by analyzing the be- 
havior of the matrices ^ (e) [see eq. (|4.26[) ] . Of special 
interest and importance is the asymptotic behavior of 
$ /. (e) at e — > which determines the long time asymp- 
totic behavior of the process. 

To demonstrate the specific features of the small- e be- 
havior of $ f. (e) we will consider two examples of con- 
trolling processes in {z}-space: exponential (Poissonian) 
and anomalous long-time-tailed. Both these examples 
can properly be treated within the diffusion model. 

a. Poissonian z(t) controlling processes. The 
Poissonian-like controlling processes can be described by 
the model of diffusive motion of the Brownian particle 
confined within the well in {z}-space. In this model at 
small e, which correspond to times longer than the time 
r z of diffusive relaxation within the well U(z) in {z}- 
space, the behavior of $/i(e) is determined by free diffu- 
sion in {y}-space and therefore 

$ f .(e) ~ e"^ 2 and w fj (t) ~ l/i 1+ "«/ 2 , (4.31) 

(j = i,t). This result shows that originally nor- 
mal CTRW processes become anomalous as a result of 
anomalous ^(^-fluctuations and the resultant behavior 
of Wj(t) = Wf j (t) coincides with that of y(i)-fluctuations, 
i.e. anomalous y(i)-fluctuations strongly modify the ki- 
netics of normal CTRW processes. 

b. Long-time-tailed z(t) controlling processes. The 
long-time tailed case is represented by the free diffusion 
model in {z}-space [U(z) = 0]. In this model the small 
e-behavior of $ f i (e) is determined by free diffusion in the 
total {y®z}-space of dimensionality n y +n z . This means 
that in the model of localized jump rates 

$ /3 .(e)~e n ^/ 2 and w f] {t) ^ l/t 1+n ^/ 2 , (4.32) 

where n yz = n y + n z . Formula (|4.32|) demonstrates 
the strong effect of (^-fluctuations on the kinetics of 
the originally anomalous CTRW processes. Noteworthy 
is that Wj (i)-fluctuations result in the increase of the 
anomaly parameter a which determines the long time 
behavior of Wj{t): Wj(t) ~ l/t 1+a . 

It is interesting to note that the effect of fluctuations 
can lead to the crucial change of fluctuation jump statis- 
tics. The fact is that usually CTRW theories assume that 
J °° dt Wp (t) = 1. This relation ensures conservation of 
normalization (or population) in CTRW processes. In 
the free diffusion approximation for the controlling pro- 
cess this conservation relation is fulfilled if the dimen- 
sionality of the space is small, when the statistics of re- 
occurrences in the jump state is recurrent. That is why 



we have assumed n y < 2 and n z < 2. These two inequal- 
ities, however, do not warrantee the same inequality for 
n yz = n y + n z . In principle, one can get n yz > 2 and 
in this case the the statistics of reoccurrences becomes 
transient which means that J™dtWj(t < 1. Such CTRW 
processes with decay are known in the probability theory 
although are not applied widely (27j . 



V. DISCUSSION 

The results obtained in this work have demonstrated 
that the Markovian representation is very useful for the 
analysis and extensions of the CTRW approach. In this 
section we will discuss some specific features of the pro- 
posed method and obtained results. 

But first we would like to emphasize the important 
point concerning the relation of this representation to 
conventional approaches applied in the theory of CTRW 
processes. The fact is that the Markovian representation 
can be considered as a convenient method of realizing 
subordination (in a fairly general form) which is conven- 
tionally used as a basis for formulation of the CTRW 
approach [T3, H3, HH. Within this representation the 
subordination is associated with the controlling process 
in {z}-space (Sec. III.B). The important advantage of 
the proposed realization consists in simplification of the 
description of non-Markovian kinetic problems by reduc- 
ing the treatment to manipulations with linear opera- 
tors. Moreover, in many cases the representation allows 
for deep understanding and modeling of real stochastic 
processes in clear physical terms. 



A. Applicability of CTRW approach 

The Markovian equation (|3.12| , applied to deriving the 
Markovian representation, is of course more general than 
the CTRW approach and therefore is quite suitable for 
the analysis of applicability conditions of this approach. 
Here we will present some comments on this point. 

The problem reduces to analyzing the statistics of fluc- 
tuation jumps controlled by the stochastic process in 
{z}-space which is determined by the evolution opera- 
tor G(e) — (e + k z + A) -1 , where k z is the jump-rate 
matrix diagonal in {x}- and {z}-spaces [see eq. (13.161) ]. 

In the limit of high localization of jump-rate matrix 
k z (in {z}-space) defined by eq. (|3.14p the Markovian 

equation (j3 . 1 2[) leads to the evolution operator G(e) cor- 
responding to the CTRW approach. In the case of ex- 
tended k z , however, the CTRW approach is not valid, 
strictly speaking. 

To find the conditions which can ensure the applica- 
bility of the CTRW approach we will consider above- 
proposed simple and fairly flexible diffusion model for 
the controlling operator A [eq. p. lip ]. 
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1) Poissonian-like statistics. Within the diffusion 
model the Poissonian-like jump statistics is described by 
the evolution operator G(e) with the Smoluchowski oper- 
ator A, in which the potential u(z) is of type of infinitely 
deep potential well, say, of width a u , i.e. u{a u ) ~ 1. 
At long times the proposed model predicts exponen- 
tially decreasing function w n (t) with the character time 
Tn ~ al/D z . 

It is easily seen that in this model the limit of high 
localization is realized for « a u , where is the 
characteristic width of the function kf . 

2) Anomalous long-time tailed statistics. The more in- 
teresting case of long-time tailed jump statistics is repre- 
sented by the free diffusion variant of the diffusion model 
(u(z) = 0), which in the limit of highly localized kf pre- 
dicts the PDF Wj(t) ~ l/t 1+ " 2 / 2 , where n z is the dimen- 
sionality of {z}-space. 

In the case of long-time tailed statistics the condition of 
applicability of the CTRW approach (applicability of the 
approximation of highly localized jump rates) is less triv- 
ial than that formulated above for Poissonian-like statis- 
tics. The fact is that in the absence of the potential one 
needs to compare the characteristic width of k z in {z}- 
space with that of the PDF <p(z, t) of the controlling pro- 
cess [see eq. (I3.9[) ] which yield the condition a\ « D z t. 

It is important to note that to satisfy this condition the 
function kf does not need to be very short range. It can 
be shown that, for example, in the case fcf ~ l/z m the 
proper parameter a& can be introduced for m > 3 [29| . 
In addition in the case m > 3 the long time behavior of 
the PDF is shown to be identical to that for really short 
range k% with the corresponding size at ■ This means that 
for m > 3 at long times « D z t the CTRW approach 
is quite applicable. 

This brief analysis shows that in the large class of 
anomalous non-Markovian models the long time behav- 
ior of the evolution operator G(t) is correctly described 
by the CTRW approach. 



In particular, let us discuss anomalous relaxation in 
highly disordered structures. In such structures the con- 
trolling processes are often quite adequately described 
by the anomalous model considered in Sec. IV.B in 
which the controlling evolution functions </>n(e) ~ £ aq 
with otq < 1 (q = 1, . . . , N) [see eq. 6H] 0. In this 
model formula (|4.15p predicts interesting behavior of the 
total controlling evolution function $(e) ~ e a °, where 

ao = YiiLi a i- I n principle, the product in equation for 
ao is convergent for properly behaving ai as a function 
of i. Corresponding criteria are known [30l ] but it is ev- 
ident that the necessary condition is a^rx, — > 1. This 
condition has fairly clear physical interpretation: if we 
escribe the controlling functions with larger numbers i 
to the structures of smaller size, then approaching of ctj 
to unity results from the evident fact that in structures 
of smaller size the anomalous effects, caused by disor- 
der of medium, are expected to be weaker. Weakness 
of effects manifests itself in the reduction of processes to 
Markovian with the increase of i, i.e. just in the relation 
cti^oo — * 1 or ai — 1 at i larger than some characteristic 
number N. 

It is also important to note another interesting predic- 
tion of the expression (|4.15[) . It shows that in the case of 
cascaded controlling process the anomaly of the processes 
in the cascade is accumulated. This effect manifests itself 
in the decrease of ao with the increase of the number N 
of coupled processes. In particular, for a large number N 
of weakly anomalous cascaded controlling processes with 
6i = 1 — 04 <C 1 one gets the value 

a^e-^^e-Jo^a-;), (5.1) 

which can be fairly small, corresponding to strongly 
anomalous process. Formula (|5.ip can shed light on the 
mechanism of formation of relaxation anomaly in disor- 
dered systems. 

2. Fluctuations of waiting time PDF matrices 



B. Extended CTRW approaches 

In Sec. IV we have studied most general features of two 
extensions: CTRW with cascaded controlling processes 
and CTRW processes with fluctuating PDFs. Below we 
will discuss in detail some particular predictions of the 
extended CTRW-approaches. 



1. Cascaded controlling processes 

The model of cascaded control is in reality a good tool 
for the analysis of the kinetic process in fractal structures. 
The cascade of coupled processes can properly model the 
kinetic coupling of structures of different size. The model 
is very useful for the analysis of kinetics of processes in 
selfsimilar and complex structures [HI]. 



The results obtained in Sec. IV.B demonstrate that 
fluctuations of waiting time PDF matrices can strongly 
manifest themselves in the kinetics of CTRW-like pro- 
cesses. It is, however, still worth to point out some spe- 
cific features of this effect. 

1) In our consideration in Sec. IV.B these fluctuations 
are assumed to result from those of the jump rate ko, 
i.e. from the dependence fco(y). For simplicity, we have 

used the model of highly localized fco(y) = £ |yt)(yt|. 
However, in accordance with conditions of applicability of 
CTRW approaches discussed above in Sec. V. A, the main 
conclusions on the effect of fluctuation in the kinetics 
of the process remain valid for delocalized dependencies 
fco(y) as well, if these dependencies are sharp enough 
(for details see Sec. IV.B). In other words the obtained 
results are valid for the wide class of models based on the 
Markovian representation. 
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2) For the sake of simplicity of the analysis, the op- 
erators A z and A y have been suggested to be indepen- 
dent of y and z variables, respectively. Nevertheless, the 
major part of conclusions of Sec. IV. B are valid in the 
case of coupled processes in {y}- and {z}-spaces (i.e. for 
[A z , Ay] 0) . The most general formulation in this case 
consists in replacement of the sum A z + A y by some op- 
erator A yz in the combined space {y <g> z). 

3) To illustrate possible modifications of the fluctu- 
ation mechanism within the above general formulation 
we will briefly discuss the simple model which describes 
fluctuations occurring at moments of jump transitions. 
In principle, this model corresponds to the highly local- 
ized z-dependence of the operator A„(z): A y (z) ~ P", 
but in the end the dependence reduces to the following 
modification of the operator C [see eqs. p.7[) and (|3.14D ] 

C = k (y)P t z t -P x P y k (y)Pf. (5.2) 

Here P y describes sudden change of y-coordinate simul- 
taneously with the jump in {a;}-space. Of course Py sat- 
isfies the normalization condition (y e |(l — P y ) = 0, where 
(ye I = Si(yj'l * s the adjoint equilibrium vector in {y}- 
space. For example, in the simplest variant of sudden 
relaxation in {y}-space P y = |y e )(y e |. 

Substitution of the expression (|5.2p into the SLE (|3.12p 
and subsequent manipulations similar to those presented 
in Sec. Ill yield 

6(e) = <ye|5 v (e)|yi>, (5.3) 

where 

G„(e) = PJe) +P VB (c)[l- W^e)]" 1 ^). (5.4) 

The parameters in this formula are similar to those de- 
fined in eqs. (|3.15p - (|3.2ip . The only difference consists 

in additional term P y in the expression for W Vj , (j = 
i,n), and in y-dependence of the parameters resulting 
from the dependence fco(y)- According to eq. (|3.17p . 
W % = P x PyW yj (e), where ?%(e) = k {y){z t \G y (e)\zj) 
and G v {e) is obtained with fco(y)- This means that for- 
mula (|5.4p coincides with eq. (|3 . 1 7|) in which P x re- 
placed by P xy — P x P y . In other words this formula 
describes jump-like migration in {x ® y}-space though 
with jump rate ko(y) which has addition dependence 
on y. Moreover, the expression (|5.4p can also be rep- 
resented in the form (|4.5p suitable for the analysis of the 
limit of diffusion-like motion governed by the operator 
£-xy — (1 — Pxy)ko(y) describing diffusion with the coef- 
ficient D(y) ~ feo(y) m some effective potential |3l| . 

It is important to note that in this model (unlike the 
general model discussed in Sec. IV. B) W(i)-fluctuations, 
caused by stochastic motion in {?/}-space, do not lead 
to the significant and universal change of the long time 
behavior of matrices w f. (t) defined in eq. (|4.24|) [see eqs. 



(|4.3ip and (|4.32p ]. Some change of iiif (t)-behavior is, 

in principle, possible due to the effect of D{y) ~ ko(y)- 
dependence, but this change is not universal, strongly 
depending on specific features of fco(y)-behavior. 

It is also worth noting that the model, which has some- 
thing in common with the considered particular variant 
of our general model, is recently discussed in ref. [22l |. 
Restricting ourselves to brief comments we would only 
like to emphasize that, as the analysis of above simple 
variant shows, the general method developed in our work 
enables one to represent the results obtained in this paper 
in very compact and general form. 

4) From very beginning in our consideration the ex- 
tensions have been discussed within the Markovian rep- 
resentation reducing the problem to the analysis of mul- 
tidimensional SLE. In so doing we have not interpreted 
the results in terms of the conventional approach ap- 
pealing to stochastic properties of fluctuating PDFs of 
consecutive jumps Wp(i) (see Sec. III). It is clear that 
the case of most pronounced effect of fluctuations cor- 
responds to strong long time tailed correlations of the 
PDFs and of course such expression can be found with- 
out difficulties. For example, one can consider the model 
assuming the process of stochastic change of Wp(t) lo- 
calized (in time) near the time of fluctuation jump. 
In this model the stochastic change of PDFs can be 
described by introducing additional fluctuation matrix 
W c (£) with which conventional convolution terms of 
type J* 2 dt{W(t — ii)W(ii — to) are transformed into 
J * 2 dhW(t - ii) J* 1 dTWdh - r)W(r - to). It is clear 
that in the limit of very short range dependencies W c (i) 
this model reduces to that considered above. However, 
even in the most general formulation it can be treated 
as a particular variant of the above-proposed approach 
based on the Markovian representation. 

The analysis of possible other variants of extensions 
which can be analyzed within the Markovian representa- 
tion will be presented elsewhere. 



VI. CONCLUSIONS 

This work concerns detailed discussion and applica- 
tions of the Markovian representation of non-Markovian 
CTRW-like processes and, in particular, non-Markovian 
CTRW-based SLE. In reality, however the Markovian 
representation, reducing the problem to the study of the 
multidimensional Markovian SLE, is more general than 
the CTRW approach and allows one not only to analyze 
the applicability of this approach but also develop some 
extensions. In our work we have considered two of them 
describing the effect of cascaded controlling processes and 
fluctuations of jump PDFs. The number of extensions is, 
however, very large. For example, interesting effects can 
be predicted in variants of CTRW approaches combining 
two above-mentioned extensions. Further generalizations 
and applications is a subject of subsequent publications. 
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